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Improving U2O2 conductance via polaron preemption: an ab initio study of Si 
doping 
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We report on ab initio electronic structure simulations of Li 2 02, where 1.6% of lithium atoms are substituted 
by silicon. It is demonstrated that this leads to the formation of conducting impurity states in the band 
gap of Li2C>2. We show that these states originate from the antibonding orbitals of the oxygen pairs and 
are remarkably stable against possible polaron formation (upon electron injection). Through this polaron 
preemption mechanism, the proposed compound is expected to show significantly higher electronic mobility 
than stoichiometric Li2C>2, which could have significant applications in lithium-air batteries. 
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In recent years lithium-air (Li-air) batteries have re- 
ceived considerable attention from both experimentalists 
and theoreticians, due to their exceptionally high theo- 
retical energy density*^ However, their practical imple- 
mentation has been severely stalled by a series of techni- 
cal challenges, including: poor cyclability, low electrical 
efficiency, and reduced capacity at high discharge rates & 
These obstacles though to partially arise from strong po- 
larization losses at the cathode material; which, in turn, 
results from the poor conductivity of lithium peroxide 
(Li202) - the main discharge product, forming at the 
cathode of Li-air batteries^^ 

Recent theoretical studies have demonstrated that sto- 
ichiometric Li2C>2 is a wide band gap insulator, and 
that charge transport in this compound takes place via 
polaron hopping4r— In their first-principles theoretical 
study, Kang et a& demonstrated that an electron injected 
into Li2 02 will be self-trapped through the formation of 
a small polaron. It was further shown that polaron for- 
mation arises due to the molecular nature of the Li2 02 
conduction band, whereby an injected electron is local- 
ized at an 0-0 site possessing an elongated oxygen pair. 
The calculated hopping rate of this small polaron was 
directly correlated to the extreemly low electronic mo- 
bility of Li2 02- At the same time, Hummelsh0j et a/ 5 
argued that Li vacancies, which are likely to be created 
during the discharge process, may induce delocalized hole 
states in the valence band of lithium peroxide, therefore 
improving Li2 02 charge transport via hole conduction. 
However, Ong et al& demonstrated that hole states may 
also become localized through the formation of "hole" po- 
larons. Although the migration barriers of hole polarons 
were determined to be lower than electronic polarons, the 
resulting conductivity of Li 2 02 was still far below tech- 
nologically desired values. Thus far it seems the only 
practically feasible way of improving lithium peroxide's 
conductivity, proposed up to the present moment, relies 
on the use of "contact states" : Zhao et aU demonstrated 
that growing of Li202 on a graphene layer induces the 
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hole-type conducting channels in lithium peroxide. How- 
ever, this approach is limited by the short penetration 
depth of surface states into Li 2 02 layers. Such states 
only penetrate only a few Li02 layers inside the lithium 
peroxide crystal, thus limiting conduction enhancement 
to only a thin surface layer. 

In this Letter we propose an alternative way of im- 
proving the conductivity of lithium peroxide. By means 
of ab initio calculations, we show that the substitution of 
a small fraction (1.6%) of Li atoms by Si impurities leads 
to the creation of additional conducting states in Li202- 
These states are shown to originate from the partial oc- 
cupation of oxygen antibonding orbitals by electrons, do- 
nated by impurity Si atoms. The elongated oxygen pairs, 
forming these states, are bound to Si impurities, and are 
not subject to polaron-induced deformations. This po- 
laron preemption mechanisim is predicted to significantly 
enhance the electron mobility of lithium peroxide. 

Our study was performed within density functional 
theory (DFT)j^ Exchange-correlation interactions were 
captured within the generalized gradient approximation 
(GGA)3 We utilized the SIESTA package^ which cm- 
ploys norm-conserving pseudopotentials and numerical 
atomic orbitals (NAOs). A flexible basis set consist- 
ing of doubled {s}, tripled {p x ,PyiPz\, and polarized 
{d} orbitals was determined essential to obtain a cor- 
rect description of the electronic structure and energet- 
ics of this system. In all calculations, structures were 
relaxed both with respect to atomic positions and unit 
cell size. Benchmark calculations on stoichiometric Li202 
produced lattice parameters of a = 3.19 and c = 7.71 A. 
This is in excellent agreement with previous all-electron 
calculations (3.18 and 7.73 performed using the 

same GGA functional, and also agrees well with experi- 
mental data (3.14 and 7.65 A)H 

We constructed a 4 x 4 x 2 supercell to model silicon 
impurities in Li2 02- This relatively large 256-atom su- 
percell allowed us to study a system with 2 silicon atoms 
per unit cell, and to calculate several possible Si distribu- 
tion patterns (discussed below) . We found it necessary to 
use a 4 x 4 x 4 fc-point mesh to ensure total energy con- 
vergence. As additional benchmark, we calculated the 
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FIG. 1. (Color online) The 4x4x2 supercell of L12O2, used 
in doping calculations. Li, O, and Si atoms are shown in blue, 
red, and green colors, respectively. The positions of Si atoms 
correspond to a uniform impurity distribution, used in the 
electronic structure analysis. 



formation energies of Li vacancies in Li2C>2. Using the 
bulk Li chemical potential as a reference, we obtained the 
formation energies of 3.39 and 3.52 eV for intralayer and 
interlayer Li vacancies, respectively. These results are in 
line with a previous value of 2.8 eV (obtained using the 
GGA functional^), as well as more recent results of 3.8 
and 4.1 eV (calculated using the HSE functional^) . 

We acknowledge that the impurity atom (Si) distribu- 
tion within a Li2 02 host matrix is a complex problem, 
that deserves a separate investigation. Many different 
configurations may exist, including possible dimerization, 
clusterization, and phase segregation. The invesigation 
of this complex problem, which should certainly be use- 
ful and interesting, is, however, outside of the scope of 
the present study. The purpose of our investigation is 
to study the influence of Si impurity atoms on the elec- 
tronic properties of Li2C>2, and to show how these impu- 
rities may improve the conductivity of lithium peroxide. 
Therefore, we have limited our study to several possible 
impurity distrubution patterns, which were chosen using 
the following assumptions: (i) we considered an "impu- 
rity regime" , with no direct interaction between Si atoms, 
by placing them at significant distances from each other; 
(ii) only substitutional Si doping in the Li02 layer was 
considered (intralayer doping) , as this configuration was 
found to be more energetically favorable than the inter- 
layer doping configuration (see below); (iii) keeping the 
Si-Si distance as large as possible, three patterns of Si 
distribution in the supercell were tested, including one 
uniform distribution, and two patterns where Si atoms 
order along certain directions in the Li2 02 structure. 

The 256-atom supercell used in our calculations is 



shown in Fig. [TJ Within this structure, Si can dope sub- 
stitutionally at two symmetrically distinct sites: within 
the LiC>2 layer (intralayer doping) and between LiC>2 lay- 
ers (interlayer doping). We have found the intralayer 
doping configuration to be 1.1 eV lower in total energy 
than interlayer doping. To further investigate intralayer 
doping, three different patterns of widely separated Si 
dopants were considered in our supercell. In the first dis- 
tribution pattern (shown in Fig.[T]), Si atoms were placed 
in two different non-adjacent LiC>2 layers along the diag- 
onal of the unit cell. This pattern models a near uniform 
impurity distribution, with a Si-Si separation distance of 
11 A. In the the second and third patterns (not shown in 
Fig. [U, we considered two possible orderings of Si impu- 
rities along a given direction in the Li2C>2 lattice: inside 
the LiC>2 layer and normal to the Li02 layer. 

Amongst all three structures, the uniform Si distri- 
bution (Fig. [T]) was found to be the most energetically 
favorable one. The next lowest-energy structure (with 
in-plane Si ordering) was found to be 0.23 eV higher in 
total energy. Therefore, only the uniform Si distribu- 
tion in Fig. [T] was considered in followup calculations. 
Lastly, we estimated the energetical stability of the Si- 
doped Li 2 C>2 structure, shown in FigJTJ by calculating its 
formation energy. Taking as references the chemical po- 
tentials of Li and Si in their bulk phases, we obtained 
the formation energy of the silicon substitutional defect 
to be -4.9 eV (per Si atom). This demonstrates that the 
occupation of Li intralayer sites by Si atoms should be 
energetically favorable during the discharge process. 

Next, we investigated the influence of Si-doping on the 
electronic structure of lithium peroxide. The total den- 
sity of electronic states (DOS) for the doped system is 
plotted in Fig. 12(a) (solid black line). The DOS for un- 
doped Li202 is also shown for direct comparison (dashed 
black line) . The electronic structure of pure lithium per- 
oxide has been extensively studied in literature^— and 
our results for the undoped case agree well with previ- 
ous investigations. The valence and conduction bands of 
Li202 are almost exclusivlely formed by the 2p-states of 
the O^ - oxygen pair. The valence band is composed of 
<7, 7r, and 7T* states, while the conduction band has a a* 
character. Fig. [2(a) predicts the following main changes 
in lithium peroxide DOS upon Si doping: (i) a new group 
of localized states is created at -7 eV near the bottom of 
the valence band; (ii) additional states appear in the gap 
between ir and ir* states in the middle of the valence 
band; (iii) new conducting states appear in the band gap 
of the system right above the valence band maximum. 

To understand these electronic structure changes, let 
us begin by investigating the structural changes induced 
by Si substitution. It can be seen in Fig. [1] that Si 
enters into a bonding interaction with 3 nearby oxy- 
gen pairs. This interaction results in the elongation of 
each O2 pair and the formation of a SiOg cluster. In- 
terestingly, the 0-0 distances in this cluster are 2.11 A, 
which is significantly larger than the inter-oxygen dis- 
tance of unperturbed O2 ions in the system (1.55 A), 
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FIG. 2. Panel (a) shows density of electronic states (DOS) 
for Si-doped L12O2 256-atom supercell. DOS for undoped su- 
percell is also shown for comparison. Panel (b) shows DOS, 
projected on the atomic orbitals of silicon and on the neigh- 
boring three elongated O-O pairs. 



and is nearly identical to the 0-0 distance of the pola- 
ronic site (2.1 A) which is formed when extra electron is 
injected into Li202^ To trace the origin of this similar- 
ity we calculated the density of states of the SiOg dopant 
cluster (see Fig. [5b). Two immediate observations can 
be extracted from this plot: (i) the states at —7 eV are 
mostly created through the hybridization of Si 3s-orbitals 
with 2p-orbitals on surrounding oxygen atoms; (ii) the 
bottom and the middle of the valence band is formed by 
the a, 7T, and ir* orbitals of the stretched O2 molecules, 
hybridized with the 3p-states of Si. However, the most 
dramatic changes take place near the band gap of the 
system. We see that the a* states of the oxygen pairs, 
surrounding a Si atom, no longer contribute to the con- 
duction band of the system. The spin-up a* states be- 
come populated and are contributing now to the valence 
band maximum of the system, while the empty spin-down 
a* states take a resonant position just above the valence 
band maximum in the band gap of the system. The sit- 
uation is strikingly similar to the one of the polaron for- 
mation in Li202, where the extra electron populates the 
tr-antibonding state of the oxygen pair, leading to its sig- 
nificant elongationji In our case the situation is very sim- 
ilar, but with an important exception: no extra charge is 
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FIG. 3. The difference between self-consistent charge density 
and the superposition of atomic charge densities, calculated in 
the plane, crossing one Si atom and one neighboring O-pair. 



The charge density is in e/A 



supplied to the system. 

Yet, the question remains: in the absence of electron 
injection, what is the source of the charge in the a* or- 
bitals? To delve further, we performed an analysis of 
charge transfer in Li202 upon Si substitution as shown 
in Fig. [3] In this plot we present the difference between 
self-consistent charge density and the superposition of 
initial atomic charges, calculated in the plane, crossing 
one silicon atom, and a neighboring oxygen pair. It can 
clearly be seen that charge is depleted on the Si site, 
and significant charge accumulation occurs on oxygen a* 
antibonding orbitals, which is consistent with projected 
density of states analysis, presented in Fig. [2b. Over- 
all, each elongated oxygen pair in an SiOg cluster has 
one electron more than a regular lithium peroxide 0\~ 
pair. This additional charge is transfered from the Si 
atom, which serves in this system as "additional source" 
of electrons. 

The most valuable effect of Si doping is, in our opinion, 
the induction of conducting states right above the valence 
band of the system. These states are formed by the a* 
antibonding orbitals of the oxygen pairs, just as pola- 
ronic states are, and essentially preempt the formation 
of polarons. Meaning they should be stable upon elec- 
tron injection, since the further deformation of oxygen 
pair is energetically quite costly. This is in contrast to 
pure Li202, where these states are located much higher in 
energy (above the band gap) and drop to form a polaron 
when an extra electron is injected. 

To verify this assumption, we performed 256-atom su- 
percell calculations with one extra electron added to the 
system. Fig. 0^ shows the DOS for the charged system. 
It can be seen that the conducting a* states are now 
partially occupied, and the DOS profile is only slightly 
modified, as compared to the undoped structure. We also 
computed the real-space distribution of electrons at the 
Fermi level of the doped system. Fig. 0b shows that the 
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FIG. 4. Panel (a) shows density of electronic states (DOS) 
for Si-doped Li2 02 256-atom supercell, doped with one extra 
electron. Panel (b) shows a spacial charge distribution at the 
Fermi level for the system with one extra electron. 



extra electron is shared between all oxygen pairs around 
Si impurities. Lastly, to examine possible polaron forma- 
tion somewhere outside the region of a Si impurity, we 
imposed a stretched O^ - pair at the maximum possible 
distance from our impurity sites. Upon structural relax- 
ation, the charged system always prefered to return to 
initial state shown in Fig. |4p. 

To complete our study we provide a first order esti- 
mate of the electron mobility fi e enhancement, that might 
be expected when electrons are injected into this doped 
system. Mobility can be calculated using the Einstein 
relation 13 fi e — eD/ksT, where ks is Boltzmann's con- 
stant, and T is a temperature. To first order the diffu- 
sion coefficient D can estimated as D = L 2 /t, where L 
is the distance between a* orbitals near neighboring Si 
atoms, and r is the electronic state lifetime. The inverse 
lifetime of electrons in a given cr*-state can be estimated 
via: 1/r = where 7 is the width of the a* DOS peak, 
and H is Plank's constant^ Substituting in the results 
from our current study, with a Si concentration of 1.6% 
(L ~ 7 A, 7 = 0.05 eV), we obtained the room tempera- 
ture mobility value of ~ 16 cm 2 /[Vs]. This value is sig- 
nificantly higher than polaronic hopping estimates calcu- 



lated for pure lithium peroxide (10 -10 — 10~ 9 cm 2 /[Vs]). 4 
The reason for this dramatic increase lies in the differ- 
ent charge transport mechanisms dominating the two 
systems: in stoichiometric Li202 charge carriers move 
through polaronic hopping, while in the proposed Si- 
doped structure charge transport likely occurs mainly 
through quantum mechanical tunneling between SiC>6 o*- 
states. However, carrier transport assisted by these Si 
dopants deserves a careful detailed study on its own. We 
speculate that even a small increase in the Si concentra- 
tion should lead to further increase of the spacial density 
of conducting <j*-orbitals, leading to the formation of an 
intrinsic conducting network inside the LI2O2 crystal. 

In conclusion, we have presented a theoretical study of 
the electronic properties of Li2 02 crystal, where 1.6% of 
lithium atoms were substituted by silicon. It has been 
demonstrated that under this type of doping additional 
conducting states are created in the band gap of the host 
crystal right above the top of the valence band. We have 
further demonstrated, that these conducting states origi- 
nate from the antibonding states of the oxygen molecules, 
which have been partially populated by charge transfer 
from the silicon atoms. These results show that this 
charge transfer elongates the oxygen pairs, and makes 
them stable against polaron formation when an extra 
electron is injected in the system (effectively preempt- 
ing polaron formation). The electron mobility of the 
proposed system is predicted to be significantly higher 
than that of intrinsic lithium peroxide, due to the non- 
polaronic nature of Si assisted carrier transport. Pro- 
vided that each Si atom creates three elongated oxygen 
pairs, we believe that electronic transport the Li2 02 lat- 
tice should be significantly improved at reasonable con- 
centrations of Si impurities. 
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